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Abstract 

The correlation spectrum of fully developed one-dimensional mappings are 
studied near and at a weakly intermittent situation. Using a suitable infi- 
nite matrix representation, the eigenvalue equation of the Frobenius-Perron 
operator is approximately reduced to the radial Schrodinger equation of the 
hydrogen atom. Corrections are calculated by quantum mechanical perturba- 
tion theory. Analytical expressions for the spectral properties and correlation 
functions are derived and checked numerically. Compared to our previous 
works, the accuracy of the present results is significantly higher owing to the 

controlled and systematic approximation scheme. 
PACS numbers: 05.45. +h 
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I. INTRODUCTION 



Typical intermittent behavior involves an irregular switching between regular and chaotic 
motion. As a result, correlation decay in intermittent systems is intermediate between that 
of regular and 'purely' chaotic systems: it follows typically a power law IQ. Hereafter 
we consider intermittency in case of one- dimensional noninvertible mappings. Correlation 
functions may then be expressed in terms of the spectral properties of the Frobenius- Perron 
operator defined by 



Here ^{x) is an arbitrary function, f{x) stands for the mapping function and a labels the 
inverse branches. The power law decay is associated with an accumulation point at the 
upper edge of the spectrum that causes serious difficulties at direct numerical studies of the 
spectral properties. Moreover, we have found that the corresponding eigenvectors tend to 
become approximately parallel, that greatly reduces the number of eigenstates which can 
be extracted from a finite size matrix representation 0. Therefore, a suitable analytical 
treatment is needed which does not involve any truncations. Such a method has already 
been applied in Refs. [@] , . Here we present another approach which is based on a suitable 
infinite matrix representation. This method has the advantage that it approximately reduces 
the problem to the well known radial Schrodinger equation for s-states. Mathematically 
significant is the fact that the original problem is transformed into an equivalent one which 
admits a nearly hermitian representation. Therefore, beyond the interesting analogy, our 
method renders possible the application of standard perturbation theory as a tool for a 
systematic calculation of corrections. This results in a rather high accuracy for the spectral 
properties, as we shall demonstrate. 

We develop our method for the treatment of the family of piecewise parabolic maps |Q 




(1) 



E 



a 




r G [0, 1] . 



(2) 
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Note that the parameter value r = 1 corresponds to a weakly intermittent situation - ) 
|T|]. It turns out that the main features of the spectrum and the eigenfunctions depend only 
on the behavior of the map near the unstable fixed point when r is close to 1, thus our 
results have a broader relevance, in fact, they bear in a sense a universal character 0. 

The paper is organized as follows. In Section II the intermittent situation (r = 1) 
is discussed, especially, contact is made with the quantum mechanical s-scattering on a 
Coulomb potential at zero energy. In Section III the nearly intermittent (r ^ 1) situation 
is considered. The problem is reduced approximately to the bound s-states in a Coulomb 
potential. Correction terms are calculated by standard quantum mechanical perturbation 
theory. Approaching intermittency, the number of the eigenstates essentially involved in 
the representation of a correlation function tends to the infinity, therefore the calculation of 
correlation functions is highly nontrivial, even if the spectral properties are already known 
0. In Section IV we present and discuss an analytic expression for correlation functions, 
that is relevant near intermittency. The detailed derivation can be found in Appendix B. 
A discussion of the results is given in the concluding Section V. Appendix A contains an 
outline and discussion of the numerical method for calculating corrections in the intermittent 
situation. 



II. THE INTERMITTENT SITUATION 

One conclusion of our numerical study done on the family of maps @ has been that 
finite size matrix approximations are poor near and at the intermittent situation, hence if 
one wants to get the spectral properties in those cases reliably, some other method is needed 
which does not involve any truncations of the matrix representations. As a first step we 
have chosen for an analytical study the intermittent map 

x' = 1 - ^\l-2x\ , (3) 

corresponding to r = 1 in Eq. (^. Using the basis (1 — (where n = 0,1,2,...), 

one obtains for the (j, A;)-th matrix element of the Frobenius-Perron-operator (|l]) the simple 

3 



closed expression 



H, 



Ak + 

V 2i y 



(4) 



These matrix elements are displayed in Fig. 1. Our aim is to find an asymptotical solution to 
the eigenvalue equation, i.e. an expression of the eigenvector for large j-values, as the failure 
of the numerical calculation implies that this numerically inaccessible part of the eigenvector 
plays an inevitable role. The matrix representation also supports this expectation as the 
largest matrix elements lie at the diagonal and decay with j as while the strip of the 
non- negligible matrix elements along the diagonal has a width of ~ a/J. Thus any truncation 
leads to a dramatic effect. Using Stirling's formula for the factorials arising in Eq. (|) for 
large j and \k — j\ ^ ^ the Hj]^ matrix element is approximately given by 

Pi I 2m2\ / /3m rrv'X / 5 1 27 16 m^\\ , , 

(where m = k — j), this expression being accurate up to the order 4. The quality of the 
approximation is demonstrated in Fig. 2, where the difference between the asymptotical 
expression (||) and the exact values @ is shown. Retaining only the Gaussian, we find the 
approximate eigenvectors 

al'^ = sin (cVk) , (6) 

and 

= cos (cVk) , (7) 

both belonging to the eigenvalue 

A^exp(-I) . (8) 

This can be proven by inserting the above asymptotic expressions into the eigenvalue equa- 
tion, replacing the summation over k with integration and evaluating it by using the saddle 
point method. The reason why we get this twofold degeneracy can be understood as a result 



of the asymptotical method which formally allows j and k to take on negative values as 
well, thus the situation is analogous with the continuous spectrum of a quantum particle 
performing unbounded one-dimensional motion. Actually these solutions get corrections 
near the origin (i.e., ai k = 0), where the asymptotic expansion of Hj^k does not hold, and, 
on the other hand they combine in a particular way to cancel each other for k < 0. As a re- 
sult, the spectrum will not be degenerated, but is continuous. Using a quantum mechanical 
analogy again, the situation resembles to a one dimensional scattering process of a particle 
on a potential containing a hard core. It is interesting to note, that the asymptotic regime 
corresponds to the laminar motion in the original x-representation, as for large k the basis 
functions (1 — x)'^'^"^^ are sharply peaked at the origin (which is now a marginally unstable 
fixed point) and nearly vanish elsewhere. The nonasymptotic regime, on the other hand, 
corresponds to the chaotic motion. Hence in this representation we are given a picture about 
intermittency that relates it to a scattering process, where the motion near the scattering 
potential corresponds to chaos. As we shall see later, in the nonintermittent situation one 
has to do with bound states in such a potential. This is in accordance with the fact that we 
have a discrete spectrum in that case. 

Nevertheless, the above asymptotic solutions are not very precise numerically. To im- 
prove this approximation, one may seek for corrections proportional to and to ^. In 
order to get them, one has to take into account corrections to the saddle point method, 
which means that one uses Eq. (^) and retains terms like those proportional to which 
give a contribution when averaging them with the Gaussian. It turns out, however, that 
these corrections are proportional to ^ and respectively, for C << 1, i.e., the correction 
terms diverge as the eigenvalue approaches 1 (which is the most interesting situation for 
us). On the other hand, the 'improved' solution shows a remarkable feature, namely, that 
for small C it depends on C and on k only through the combination C'^k: 

a? = KC'k) . (9) 

If one demands this dependence from the beginning and assumes that C<<1, k » 1 



while no assumption is made about C^fc, then one gets up to order i in j and to order 
in C (defined now through Eq. (|D) the equation 

h'\z) + l-h{z) = , (10) 

{z standing for C^j) whose two hnearly independent solutions are (cf. [p!0| ) 

h^'\z) = ^J,{^) (11) 

and 

h^'\z) = ^Y, , (12) 

Ji and Yi standing for the first order Bessel and Neumann functions, respectively. Eq. 
(PUP is identical with the radial Schrodinger equation for s-wave scattering in a Coulomb 
potential at zero energy. As discussed above, it is valid far from the 'chaotic core' of the 
effective scattering potential which for small values of k gives rise to corrections. We shall 
see, however, that the range of validity of Eq. ([ToD extends down to k = 1 when C 0. 
The boundary conditions follow from the restrictions that the eigenfunctions in the original 
'x-space' should have an integrable singularity at x = and their integral over the whole 
[0, 1] interval should vanish. This latter follows from the orthogonality of the left and right 
eigenfunctions (belonging to different eigenvalues) and from the fact that the identically 1 
function is the left eigenfunction of the Frobenius- Perron operator for C = (i.e., for A = 1). 
Note that it is valid only in the case of permanent chaos and does not apply for transient 
chaos. The eigenfunction of the Frobenius-Perron operator is given by 

oo 

M^)=j:9A^'^y'^' , (13) 

j=0 

where the asymptotical form of gj is given by some linear combination of h^^^ and h^"^^ (see 
Eqs.([Tl|),(p]^)). As the asymptotical form of both the Bessel and the Neumann function is a 
phase-shifted sinus with a one-over-square-root type amplitude, for fixed C and increasing j 
the moduli of the coefficients gj grow like j^. This growing is, however, superimposed by the 

6 



exponential decay of the (1 — x)"^^^^ factor, making the infinite sum in Eq. ([I3D absolutely 
convergent for < x < 1. One can even show that the eigenfunction (pd^) is integrable 
near zero (one estimates the summation with an integral, after a term-by-term integration 
over x). This is true essentially because Y.'jLoQj I {'^3 + 2) is finite. One can express the 
requirement of the integrability of the eigenfunction near zero as a boundary condition to 
Eq. ([lO|) stating that /^°° h{z)/z is finite (the lower limit of the integration range being an 
arbitrary positive value). Thus we can say that both propagating solutions ( |TT1) and ( |T2|) of 
Eq. ([ToD are allowed by the boundary condition at the infinity, as is usually the case in the 
customary scattering problems. 

The next issue is how one can determine the actual eigenvector gj for a given (small) C 
knowing the asymptotical solutions h^^\C'^j). The asymptotical solutions are actually rather 
accurate approximations even for small j-s when C is small. Nevertheless, at or near j = 
they do not satisfy the eigenvalue equation, thus giving rise to a correction term. If one starts 
with the proper linear combination of the two asymptotic solutions, then the correction term 
rapidly decays. This is the condition which selects the proper asymptotics. (If, however, not 
the proper linear combination has been chosen, then the correction term itself contains an 
asymptotic part, describing a 'reflection' at the origin j = 0.) The numerical procedure is 
outlined and discussed in Appendix A. According to those considerations, for small C, i.e., 
for an eigenvalue near unity the eigenfunction has the form 

Mx) = 2{l-x)-2f^cfjjJc^^){l-x)'^+' . (14) 

i=o ^ ^ 

(Here we have multiplied by —2 in order to exhibit the similarity to the C = case when 
just the first term of the r.h.s. of Eq. (p^ ) remains.) The sum may be evaluated by replacing 
it with an integral. Neglecting terms proportional to times a nonsingular function (which 
come from the derivatives of the summand at j = when the Euler-Maclaurin Summation 
Formula is apphed) one obtains (cf. |T^) 
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For C — s> the second term vanishes everywhere except near x = 0, where it has a sharp peak 
with unit area. As the first term is nothing but the normahzed stationary probabihty density 
of the map (which is at the same time the eigenfunction of the Frobenius-Perron operator 
with unit eigenvalue (i.e., C = 0), one can see in which sense the hmiting eigenfunction is 
approached and also, that the integral of the eigenfunction over x vanishes. 

One can see now that the spectral properties around the upper edge of the spectrum 
depend predominantly on the laminar motion of the map near its marginally unstable fixed 
point, as the influence of the other parts of the map can show up itself through correction 
terms, which, however, become negligible when the eigenvalue approaches unity. On the 
other hand, corrections for larger values of C can be calculated with relatively little numerical 
efforts as the nondecaying asymptotical part of the eigenvectors (which has been previously 
the root of the numerical difficulties) is already taken into account analytically. 

III. THE ASYMPTOTICAL SOLUTION OF THE FROBENIUS-PERRON 
EIGENVALUE EQUATION NEAR INTERMITTENCY 

We seek for the asymptotical form of the matrix elements Hj j^ of the Frobenius-Perron 
operator H corresponding to the family of maps (§) on the basis 

Cn{x) = {fr')'{x) (l-2/r^(a;))'" 
= (^^^-ra;) (l-x)'"(l-rx)'" (n = 0, 1, 2, ...) (16) 

This form is suggested by the symmetry of the map as well as by the structure of the 
Frobenius-Perron operator (p. Indeed, for a symmetric fully developed chaotic map defined 
on the interval [0, 1] the Frobenius-Perron operator can be expressed as 

•1 n „ w .\\ n n 



(h^) (x) = {fr')' (x) (2 - (2 - fr"(^))) + (2 + (2 " ^'"'^"^0) 



(17) 



which implies that any polynomial v^(x) goes over under the application of H into a lin- 
ear combination of the functions (/r^) (^) (l ~ 2/^^^(a;)^ which, if they are themselves 



polynomials, constitute a suitable basis. In our concrete example they coincide with those 
given above (Eq. (p!6|)). One can see that this basis goes over for r ^ 1 to that used in the 
intermittent situation. 

Introducing the functions 

K (x) = (HCk) (x) , (18) 
the matrix elements Hj k are defined by 

oo 

hk {x) = E H,,kQ ■ (19) 

As one has to do here with polynomials, a numerical evaluation of the matrix elements for 
any given indices j, k is not difficult, however, unlike in the intermittent situation, a closed 
analytical expression (a counterpart of Eq. (^) this time does not exist. Nonetheless an 
asymptotical expression can be derived. Our starting point now is Eq. ([T9|), where we insert 
the expressions of the basis functions (|1^). Introducing the variable z = 1 — 2fj^^{x) = 
(1 — a;)(l — r x) we get 

2-{4fc+i) ^i^^^^i^ ^)2fc + 
^2-(4fc+i) (1 _ ^ ^) (1 _ ^)2fc ^2 - r - r zf" (20) 

oo 
j=0 

The exact matrix elements may be calculated by comparing the coefficients of the polyno- 
mials in both sides. The result is shown in Fig. 3. 

In order to derive an asymptotical analytical approximation, we extend the variable z to 
the whole complex plane and determine the matrix elements by Cauchy's formula as 

Hj,k = ^Jhk{z)z-^''^'\ (21) 

where the integration contour encircles the origin and hk{z) stands for the left hand side of 
Eq. (^) . For large j and k values the integrals in Eq. (|21|) (each involving one of the terms 
of hk{z)) may be evaluated by the saddle point method. It turns out that for a given j the 
matrix elements have a maximum versus k, namely 
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TJ ^ 



2(l + 2r-r2) 2{l + 2r-r^)\ VJ 



It has been assumed that — is of order unity (i.e., the expression ( ^2]) is vahd near the 
maximum of the matrix elements in a strip of width ^ y/J)- The most important difference 
between the intermittent and nonintermittent situations is that the maximal matrix elements 
lie now above the diagonal, and, when increasing the index j, their distance from the diagonal 
increases faster (linearly), than the spread of the significant matrix elements around the 
maximum (square root type increase). As a consequence, as we shall see, the eigenvectors ak 
decay for large k exponentially. Therefore, we have to evaluate Eq. ( pT]) not along k ^ y+t-^' 
as we would expect from Eq. (|22|), but along another line k ^ bj. The quantity b as well as 
the asymptotics oc exp {—a k) of the eigenvector is to be determined from the condition 

Hj^k exp {-a k) oc exp {-a j) , (23) 

k 

where the dominant terms of the sum come from the k values near bj. Explicitly, we get 

■ ' (24) 



3 — r 
and 



a = 2 In ( - — - ] . (25) 



Evaluating the expression of the matrix element Hj^k by the saddle point method around 
k = bj in a. strip of width y/j, we get 

2(l + 2r-r2) 3-r \r J ^ \r Jx 2{l + 2r-r^)l 



X 
771 



where 



X ( 1 + p + 6i,3 p^) X + (62,0 + h,2 / + &2,4 + &2,6 /) 



^ _ (r^ - 12 + 36 - 28 r - 9) (3 - r) 

~ 4 (l + 2r-r2)2 ^ ^ 
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^ _ (r^- 16r2 + 24r + 3) (3-r)^ 

r^-33r6 + 235r5- 687r4 + 835r3-171r2- 243r-57 , , 

02 = T (29) 

48(l + 2r-r2)3 ^ ' 

(r^ + 8 - 240 r'^ + 1496 - 4010 + 4664 - 1216 - 952 r - 183) (3 - r)^ 

"2,2 ^„ , ^ ^74 (30) 



32 (l + 2r-r 

^ _ (4r^-9r6-122r^ + 629r^- 1040r3 + 459r2 + 186r + 21)(3-r)^ 

~ 24 (l + 2r-r2)^ ^^^^ 



^ _ (3-r)^(r^-16r2 + 24r + 3)^ 



and 



a; = (33) 



P=^f^. (34) 



Figs. 4, 5 show the differences between the approximate expressions (p2|) , (^) and the 
exact matrix elements. 

Our next taslc is to solve the eigenvalue equation 

oo 

Hj,kSk = . (35) 

fc=0 

We approximate the summation by an integral (considering to be a continuous variable), 
extend its lower limit to — oo, insert the asymptotical expression (^) for the matrix ele- 
ments and evaluate the left hand side by the saddle point method. In order to make it 
systematically, we write 

Sk = (p{k) exp (—a k) , (36) 

and assume that (f{k) grows (or decays) asymptotically at most like a power of k. Then 
we insert the expression (^) of Hj ^, expand (p{k) around bj on the l.h.s. up to fourth 
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order, and perform the integration. The result can be written down for arbitrary r, but 
for simphcity it will be presented for r = 1 — e, where e << 1, i.e., near the intermittent 
situation. As the l.h.s. of the equation becomes a function of bj, while at the r.h.s. ip{j) 
stands, the latter should also be expanded around bj. Introducing the new independent 
variable 



X = 4:ebj 



(37) 



and writing 



(38) 



we arrive at the equation 



eX 



+e 



0. 



(39) 



For small e the term proportional to e may be neglected in the first approximation. The 
resulting equation. 



xC - + = 

eA 



(40) 



has a solution which diverges for a; ^ oo slower than exponentially (actually, as a power) 
only if is a positive integer n. In this case 



^ — — X L\_i{x) 



(41) 



and the eigenvector Sj^^ (cf. Eqs.(^)-(|38|)) is expressed as 



sj") = Aebj Ll__-^{4ebj) exp(-aj) , 



(42) 



where Ll^_-^^{x) stands for the generalized Laguerre polynomial |IT3|. The eigenvalue A is 
given (up to first order in e) by 



A„ = 1 - -n. 



(43) 
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Let us consider the transition to the intermittent case. Eq. (^31) imphes that in that case 
(i.e., when e ^ 0) the spacing between neighbouring eigenvalues vanishes, thus we get a 
continous spectrum. In order to get the eigenvectors, we fix the value of A, thus also of en 



and then take the limit e 0, or, equivalently, n oo. We obtain |12 



4 = Ito { fc:^ 1 e-' = v/8(l - X), J, I 2^/8(1 - A), ) (44) 



n \ n 



Provided that 1 — A << 1, we may write according to Eq. (P) 1 — A ~ thus we get 

C 



'"^J\/JJ^[C^/J) , (45) 

in accordance with Eq. (|TT]), which gives the dominant contribution near the upper edge of 
the spectrum (cf. the discussion after Eq. (|I^)). 

When deriving Eqs. we have made use of the same argument as that applied 

at the solution of the radial Schrodinger equation with a Coulomb potential. It is indeed 
possible to cast Eq. (^Of ) to the same form. To do this, let us introduce the new independent 
variable 

p = Kx (46) 



and the new function 



X(p)=eexp(-^) , (47) 

where 

Then we get from Eq. ( ^0|) 

~x" --x = --r^x- (49) 

This is the well-known form of the radial Schrodinger equation. Note that the actual radial 
wave function R{p) corresponds to x/P- The term — ^ represents an attracting Coulomb 

13 



potential and E = is the Rydberg- formula. As usually, bound states are associated 

with positive integer values for k. We may apply the transformation (^6])-(0) also to Eq. 
(p9|). Before doing that we simplify somewhat the correction term. We shall be interested 
only in the next order correction to the eigenvalue (^3|), thus in the order e term of Eq. 
(p9|) we may express all the higher derivatives in terms of C,' and of ^, as we may assume 
at the given accuracy that ^ satisfies Eq. (^O]). Thus we get for the n-th eigenvalue and 
eigenfunction 



n[n 



Performing now the transformation 



, we get 



0. 



(50) 



- X 



P 



-X - e 



3n 1 nl 1 \ fnl l\ , 



;X- 



(51) 



Applying first order perturbation theory, we may express the first correction of — as the 
diagonal matrix element of the perturbing operator, i.e.. 
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oo 

-e / dpxnip) 



3n 1 nl 1 \ fnl 1\, 



(52) 



Using Eq. (^8]) this leads to the expression of the eigenvalue A 



n 



A = 1 - -e 
2 
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n n^ + 0{e 



exp 



e 7 . 
-nl - + -t 
.2 8 



0(e3 



(53) 



which is valid in e up to second order, provided that not only e, but also en is much smaller 
than 1. A comparison of these results with those of a direct numerical measurement is 
displayed in Fig. 6. 

As for a comparison in case of the eigenf unctions, first we have to determine them in 
coordinate space. Eqs. ( ^^ and ( [T^ ) imply 



j ^1^2 1 (4e6j ) exp ( - aj ) 



r + 1 



rx\ [1 — xS^Hl — rx 



(54) 



where 
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k{x) = ^[a-2ln{{l-x){l-rx))] . (55) 



Using Rodrigues' formula 



n\ X dx 

we may evaluate the integral in (|5^) to get 



L«(x) = l--^(x"+V^) (56) 



fnU ^ n fr + l \(/t(x)-l)"-i 



4e6 V 2 J k(x)"+i 

As our asymptotical method neglects any effects which appear at j = 0, the integral of these 
functions does not vanish. Therefore, a suitable multiple of P{x) = 2(q{x) must be still 
subtracted. The resulting expression, however, gives a reasonable fit only quite near to the 
intermittent case r = 1, so that for r < 0.99 one should calculate corrections as well. The 
essence of the problem is not to push perturbation theory to higher order, but to take into 
account the second independent solution of the equation (|39|) and to add corrections near 
the origin, where the asymptotic expansion does not hold any longer. The procedure is in 
complete analogy with the intermittent situation (cf. Appendix A). It may be surprising 
that the second independent solution also plays a role. Indeed, at this point the analogy with 
the quantum mechanical case breaks down. The reason is that the physical meaning of the 
eigenfunctions is different, especially, they are differently related to probability distributions. 
Therefore in the quantum mechanical case a singularity in the origin is not allowed, while it 
is allowed in our case. Note that both solutions correspond to the same eigenvalue, i.e., the 
boundary condition at infinity completely determines the spectrum. The calculation of the 
second type of eigenfunctions is simpler after a Laplace transform, which brings us back to 
the original coordinate x. Indeed, 



x) = Y.U^ebj)e--^Cjix 

j=0 



°° /r + 1 

^ en(4e6j>-"^' -rx){l- xf^{l - rxf^ (58) 



^ 2 
'r + 1 \ 1 



\ I j-oo 



■k{x)z 



2 
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where k{x) is given by Eq. (^). Taking the Laplace transform of Eq. (^) we get 

en(0) + {n + l- 2k)^{h) - k{k - = , (59) 

where 

= / dzUz)e-''' (60) 



stands for the Laplace transform of Choosing ^„(0) = in Eq. (|59|) and substitut- 

ing (^) we arrive at the 'regular' eigenf unctions (|5^) again. The 'irregular' (logarithmic) 
eigenfunctions correspond to the choice ^„(0) 7^ 0. Choosing ^„(0) = —1 we get 

Taking the inverse Laplace transform we get the second independent solution for ^n{_z\ It 
can be given as an infinite sum, but we do not reproduce here this cumbersome expression. 

In order to get a reasonable approximation for the eigenfunctions even for relatively 
small r values, we apply the scheme decribed in Appendix A (cf. Eqs. (|A1|) - (|A5|) ) . Now 
the two independent asymptotical solutions /i^-^'', correspond to the regular and irregular 
solutions for ^^(2), and the exact matrix elements Hj^h valid in the nonintermittent case (cf. 
Eq. (0)) should be used. As a demonstration, in Figs. 7-9 the results for the first three 
eigenfunctions at r = 0.95 are compared with those of the numerical measurements. In this 
case the corrections played already an important role, 37 terms had to be included. We also 
found that it was essential to use good approximations for the eigenvalues (cf. Eq. (|53D). 

IV. CORRELATION FUNCTIONS 

Recalling the definition ([l|) of the Frobenius-Perron operator ff, one may cast a correla- 
tion function 



dxB(^f^'\x)) A{x)P{x) (62) 

(/W(a;) standing for the t-th iterate of the mapping f{x)) to the form 

16 



Cf'^ = dx B{x) H' {A{x)P{x)) i . (63) 

Expanding now A{x)P{x) in terms of the eigenf unctions of the Frobenius-Perron operator 
the action of reduces to a muhiphcation of each term by the t-th power of the corre- 
sponding eigenvalue. The details of the calculation outlined here are presented in Appendix 
B. There we make use of the previously introduced infinite matrix representation ( p!9[ ) and 
apply Eq. ( ^2] ) for the eigenvectors. This is allowed under the assumption that we are close 
to an intermittent situation (i.e., corrections and the irregular Coulomb functions may be ne- 
glected) and also that the correlators A{x) and B{x) have already zero mean with respect to 
the natural measure, therefore, a substraction of a multiple of P{x) from the eigenfunctions 
is not necessary. Then we apply identities for the generalized Laguerre polynomials, substi- 
tute sums with integrals (which is made possible again by the closeness of the intermittent 
situation), and finally end up with the expression 



1 



r (4e6) 



2 



Cf'^ ^ I dxP(x)B(x) 2 (64) 

Jo ^' ^ ^ [4e6r + (l-r)(a-21n((l-x)(l-rx)))]' ^ ' 

(f (- - iexD ( ^ ~ ^^'^ - «) - 2 In ((1 - a:)(l - rx)) (Aeb - ajl - r)) \\\ 
y \2 2 2 4e6r + (1 -r) (a-21n((l -x)(l -rx))) )))' 

Here 

r = exp(-(i + I.^)<) (65) 

e = 1 - r (66) 

a = 21n(^) (67) 

b = (68) 
3 — r 

A comparison with the results of a direct numerical measurement of the correlation function 
is displayed in Fig. 10. (cf. Fig.l. in Ref. 0, where the same is displayed with a rougher 
analytical estimate). The functions A{x) and B{x) are given by 

1 if X < 5 X 10-2 
A{x) = B{x) = { (69) 

otherwise 
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V. CONCLUSION 



A new analytical method has been developed for the determination of spectral proper- 
ties near and at the intermittent situation of fully developed chaotic one-dimensional maps. 
We made use of an infinite dimensional matrix representation, determined the asymptoti- 
cal expression of the matrix elements and reduced the eigenvalue equation to a differential 
equation, which in the lowest order coincides with the radial Schrodinger equation of the 
hydrogen atom for s-states. This allows for the calculation of corrections by using the stan- 
dard quantum mechanical perturbation theory. We have demonstrated this by calculating 
the correction to the eigenvalue. The eigenfunctions were also analytically determined and 
compared with the results of numerical measurements. Finally, we derived an analytical 
expression for the correlation functions. It has been shown to fit well the numerical data. 

One may think that our study applies only to a very specific example. Actually, the 
method may be applied to any fully developed chaotic one-dimensional map which is near 
to an intermittent situation and which has analytic inverse branches. The universal character 
of our results (within the class of mappings mentioned above) has been discussed in Ref. 

We expect that the method can be extended to the case of repellers, where one also 
has a complete symbolic dynamics. As a motivation we mention that in case of repellers 
a special interplay of transient chaos and intermittency appears, as discussed in the recent 
papers [0, ||19| which the interested reader may consult. 

APPENDIX A: CALCULATING CORRECTIONS IN THE INTERMITTENT 

SITUATION 

The numerical procedure is the following: one assumes that for some n the correction 
term 5gj is negligible when j > n, and then determines both the proper asymptotics and 
6gj from the requirement that the sum 

n 

E (H)' (Ai) 

j=n/2+l 
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be minimal (the elements in the second half of the correction vector are involved). Explicitely, 
it gives 

6gj = cos (a) 5^ ^ + sin (a) (A2) 

and the proper asymptotics is given by 

cos (a) hf + sin (a) hf , (A3) 

the accurate eigenvector gj being the difference of the above two quantities. Here 6^^ is the 
solution of the equation 



fc=0 k=0 



These are the first n of the set of the eigenvalue equations. Note that if the correction 
vector 6gj is of length n (as is assumed here), its last n/2 components still enter the next n 
equations as well. Thus, when these components are small (as required), the error caused 
by them is small and the procedure is consistent. Indeed, for C = 0.1 and n = 10 these 
components are less than 10~^ (cf. Table 1). The 'phase shift' a is given by 



^(1)^(2) 

tan{2a)= --^=:^:,r;j l,^,^^ , (A5) 



^2^j=n/2+l"j "j 



as implied by Eqs.([Al|), ([A^). Numerically we get, as Table 1 demonstrates, that in case 
of C — > the 'phase shift' vanishes and 5gj ~^o,j- This result can be understood as 
follows. For C «1 and for < j << ^ the asymptotical solutions h^'^{z) can be written 
approximately as 

h^^\z)^-- , (A6) 

TT 

hence the eigenvector gj is of the form 

1 2 

gj = -C'^j cos a sin a + 6gj . (A7) 

2 71 
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As the matrix elements Hj^k are of order unity (compared to C), if a = O(C^), then Sgj will 
be also of order C^. For j > the asymptotical solutions are already accurate to order 
(provided that C is small enough), thus for those values of j 6gj practically vanishes. In 
view of Eq. ( [A4| ) this means that the above estimates for a and 6gj indeed hold. Note that 



for j = the combination (1 — X)6go ~ ^^9o enters the eigenvalue equation, hence (unlike 
when j > 0) 6 go = 0(1). Another important observation refers to the boundary conditions 
to Eq. (p!oD : as we have seen, for C « 1 we get h{z) oc z (cf. Eqs.(|lT]) and (|A3|)), just like 
in the case of the quantum mechanical scattering on a pure Coulomb potential. 



APPENDIX B: DERIVATION OF AN ANALYTIC EXPRESSION FOR THE 

CORRELATION FUNCTIONS 



introducing 



and 



a 



A,B 



dx B{x) {A{x)P{x)) 



/ dxB{x)Hnj2a{j)Q{, 
•^0 \j=0 
«1 oo oo 

dxB{x) Y.fij)E{H\Xkix) 



E 



j=0 k=0 

dx B{x) Ck{x] 



p(x) = dx B{x) C^ix) 

Jo 4ei, 



we may write Eq. (|B1D as 

Cf'^ ~ 4^ / Jo dy(3{x)Gt{x, y)a (^-^ 
Making use of the eigenfunctions Eq. (|4^) we may write 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 
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Jo 4S ^ ^fcj 4^ 



(B6) 



Let us represent a [^) (cf. Eq. (P^D) as 



a 



y ^ s OO OO 



where 



(B7) 



(B8) 



Then we may write 



dyGt{i,y)al^-^^ 



1 - ^) y ) i.li(y')AfxLt,(x) exp (-^x ) (B9) 



= /^"^rfyexp (- (1 - ^) y) xexp (-^x) £ Lt,(£)iLt,(y)A*a 



that imphes 



(BIO) 



Here A; is of the form exp(— 7/) (cf. Eq. (|53D) where 7 = | + |e^. Inserting this and applying 
the identity [|1^ 



In 



(Bll) 



we obtain 



G,{x,y) = I exp - " ^5) ^) E^^i(5^) Ki(^) " ^?(^)) 

where r stands for exp(— 7^). Let us insert the contour integral representation WE 



(B12) 



c z — y \z — y 



(B13) 
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where the contour C encircles the point z = y. We get 



G,{x,y) = -—exp -^{x-y)] i dz^^^Y L} ,{ 
*^ '^^ 2m ^\ 4eb^ ^'Jjc {z - yf fr[ '"^ 



/-I 



Z T 



\z - yj 

Using the generating function of the Laguerre polynomials [l7|, i.e., the identity 

1 



(B14) 



1=1 



that holds true for 



(i-uy 



\u\ < 1 



exp 



xu 
u-1 



(B15) 



(B16) 



we get finally 



Gtix.y) = —- exp ( rix — y)] (p dw , , 

^ y2m ^\ 4ey ^'Jjc (w(l - r) - 1)2 



exp ^-yw - 



(B17) 



Here the new complex variable w = z/y has been introduced. The contour C encircles 
accordingly the points w = 1 and w = 1/(1 — r), the latter coming from the condition (|B16D . 
Inserting Eqs.( |B17| ), (|B3|) , ( [T6| ) into Eq. ( P5| ) the integration over x may be performed to 
get 



Cf'^ ^ 2ebT dxP{x) B{x) 
Jo 



X 



dw- 



a {Aebw — a) 



2m J Jc [w [4e6r + (1 - r) {a - 21n((l - x){l - rx)))] - {a - 21n((l -x)(l - rx)))X 

(a{l - r) {Aeh - a) - 2 In ((1 - x)(l - rx)) (4e6 - a(l - r))\ 



(B18) 



8eVr / dxP(x)B(x)d 



4e6r + (l-r)(a-21n((l-x)(l-rx))) / 
X [4e6r + (1 - r) (a - 2 In ((1 - x){l - rx)^^''^ 

Here we have also used the fact that in our case P{x) = 2 (/r^) (^)- Furthermore, 



a{p) = j^ dyexp{-yp) ^ 



(B19) 



Consider now the expression 

roc 

-d'{p)= dy exp {-yp) a{y) 
Jo 



X]a(i)exp(-^)j) . 
i=o 



(B20) 
(B21) 
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Comparing this with the expansion (cf. (pi] 



(JXJ LXJ -1 

A{x)p{x) = = E^Uhn^) (i - 



X 



j=0 



j=0 



2j 



(B22) 



we get 



a'{p)^2A[fC- (l-exp (-| 



(B23) 



Inserting this into Eq. ( [B18|) we get the final expression for the correlation function: 



A,B 



1 r {AehV 

dx Fix) B(x) 

^ ^ [4e6r + (l-r)(a-21n((l-x)(l-rx)))]^ 



(B24) 



xA\f\- — - exp , 
V \2 2 ^ \ 2 



1 a(l - r) (4e6 - a) - 2 In ((1 - x)(l - rx)) (4e6 - a(l - r))^ 



4e6r + (1 -r) (a - 2 In ((1 - x)(l - rx))) 
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TABLES 



c 


1.0 


0.5 


0.1 


0.05 


A 


0.96923323 


0.99221794 


0.99968755 


0.99992188 


a 


-0.26334376 


-0.04975108 


-0.00103697 


-0.00024860 




-0.55711388 


-0.85178915 


-0.99357189 


-0.99831575 




0.09131282 


0.01196515 


0.00000560 


-0.00000357 


S92 


0.05735793 


0.00736996 


-0.00001951 


-0.00000816 


^93 


0.03704115 


0.00553255 


0.00003283 


0.00000624 


S94 


0.01757775 


0.00255459 


-0.00000197 


-0.00000159 


S95 


0.00679480 


0.00111081 


0.00000318 


0.00000036 


S9e 


-0.00025431 


-0.00000901 


-0.00000048 


-0.00000012 


^97 


-0.00374950 


-0.00063731 


-0.00000183 


-0.00000020 


S98 


-0.00472257 


-0.00087973 


-0.00000216 


-0.00000018 


S99 


-0.00406673 


-0.00084036 


-0.00000206 


-0.00000016 



TABLE I. Dependence of the phase shift and the correction coefficients on the parameter C 
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FIGURES 



FIG. 1. The exact matrix elements Hj ^ in the intermittent case. 

FIG. 2. The difference between the analytical approximation and the exact values 
for the matrix elements Hj k in the intermittent case along the diagonal. 

FIG. 3. The exact matrix elements Hj j^ in the nearly intermittent case (r = 0.7). 

FIG. 4. Solid line: the difference between the analytical approximation ( ^2|) and the 
exact values for the matrix elements Hjk in the nearly intermittent case r = 0.7 along 
the maximum k = -^j- Dashed line: the same for the difference between the analytical 
approximation (^) and the exact matrix elements. 

FIG. 5. Solid line: the difference between the analytical approximation ( p2|) and the 
exact values for the matrix elements Hj k in the nearly intermittent case r = 0.7 along the 
line k = -^j- Dashed line: the same for the difference between the analytical approximation 
(p6|) and the exact matrix elements. 

FIG. 6. Comparison of the analytical expression (|53| ) with the results of a direct nu- 
merical measurement for the first three nontrivial eigenvalues near r = 1. Diamonds: Ai, 
crosses: A2, squares: A3 (numerically measured data). The dotted and solid lines display 
the corresponding analytical results up to first and second order in e, respectively. 

FIG. 7. Comparison of the analytical expression (|57|) with the results of a direct numerical 
measurement for the first nontrivial eigenfunction s^^\x) at r = 0.99. Solid line: analytical 
approximation, dashed line: numerical result. The eigenfunction is normalized such that 
the integral of its modulus is unity. 

FIG. 8. Same for the second eigenfunction s^'^\x). 

FIG. 9. Same for the third eigenfunction s^^\x). 

FIG. 10. The ratio between the numerically determined correlation function and the 
analytical expression ( |B24| ) at three control parameter values. Solid line: r = 0.9999, 
dashed line: r = 0.99, dotted line: r = 0.98. The functions A{x) and B{x) are given by Eq. 
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